Binding of ATP to the CBS domains in the C-terminal region of CLC-1 
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The common gating of CLC-1 has been shown to be inhibited by intracellular adenosine triphosphate (ATP) in 
acidic pH conditions. Such modulation is thought to be mediated by direct binding of ATP to the cystathionine 
p-synthase (CBS) domains at the C-terminal cytoplasmic region of CLC-1. Guided by the crystal structure of the 
C-terminal domain of CLC-5, we constructed a homology model of CLC-1 's C terminus and mutated critical amino 
acid residues lining the potential ATP-binding site. The CLC-1 mutations V634A and E865A completely abolished 
the ATP inhibition of CLC-1, consistent with the loss of ATP binding seen with the corresponding mutations in CLC-5. 
Mutating two other residues, V613 and V860, also disrupted the ATP modulation of CLC-1. However, placing aro- 
matic amino acids at position 634 increases the apparent ATP affinity. Mutant cycle analyses showed that the modu- 
lation effects of ATP and cytidine triphosphate on wild-type CLC-1 and the V634F mutant were nonadditive, 
suggesting that the side chain of amino acid at position 634 interacts with the base moiety of the nucleotide. The 
mutation effects of V634F and V613A on the ATP modulation were also nonadditive, which is consistent with the as- 
sertion suggested from the homology model that these two residues may both interact with the bound nucleotide. 
These results provide evidence for a direct ATP binding for modulating the function of CLC-1 and suggest an overall 
conserved architecture of the ATP-binding sites in CLC-1 and CLC-5. This study also demonstrates that CLC-1 is a 
convenient experimental model for studying the interaction of nucleotides/ nucleosides with the CBS domain. 



INTRODUCTION 

The resting conductance of the sarcolemma of skeletal 
muscles consists of a CI component largely contrib- 
uted by the CLC-1 Cl~ channel (Jentsch et al., 2002, 
2005; Chen, 2005). It has been shown that the inhibi- 
tion of the sarcolemmal CP conductance is critical for 
reducing the shunting conductance and, therefore, for 
helping the muscle overcome fatigue resulting from the 
accumulation of extracellular K + ions after continuous 
firing of action potentials (Pedersen et al., 2004, 2005). 
This antifatigue mechanism is thought to result from 
the ATP inhibition of the common gating of CLC-1 in 
the presence of an acidic intracellular pH (Bennetts 
et al., 2007; Tseng et al., 2007) . Presumably, ATP acts on 
the protonated CLC-1 channel by binding to the cystathio- 
nine (3-synthase (CBS) domain located at the C-terminal 
cytoplasmic region (Meyer and Dutzler, 2006; Markovic 
and Dutzler, 2007; Meyer et al., 2007) , although the mo- 
lecular mechanism of the CLC-1 modulation by ATP ap- 
pears to be unsettled. Although our study (Tseng et al., 
2007) and that of Bennetts et al. (2005, 2007) demon- 
strated an ATP inhibition of the common gating of 
CLC-1, others have suggested that ATP may exert its 
actions on CLC-1 indirectly via an accessory protein 
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(Zifarelli and Pusch, 2008) . We have recently shown that 
the controversy may result from different degrees of 
oxidation of the CLC-1 channel molecule in various 
experimental conditions (Zhang et al., 2008). However, a 
study to provide evidence for a direct binding of ATP to 
the predictive ATP-binding site of CLC-1 is imperative. 

So far, the high resolution structure of the C-terminal 
domain of CLC-1 is not available. However, the C-terminal 
structure is known in several CLC molecules, including 
CLC-0, CLC-Ka, CLC-5, and CmCLC (Meyer and Dutzler, 
2006; Markovic and Dutzler, 2007; Meyer et al., 2007; 
Feng et al., 2010). All these C-terminal structures show 
conserved CBS domains. The CLC-5 's C-terminal struc- 
ture is particularly helpful because it is the only CLC 
structure in which the ATP molecule can be seen at the 
ATP-binding site in the CBS domain. Mutating two 
CLC5's amino acid residues that interact with the 
bound ATP, as guided by the x-ray structure, reduces 
the ATP binding affinity (Meyer et al., 2007) . Based on the 
C-terminal domain of CLC-5, we constructed a homol- 
ogy model of the C terminus of CLC-1 . We mutated several 
previously unexamined amino acid residues in CLC-1 
that appear to interact with the bound ATP molecule. 
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The results show that mutating these residues greatly 
alters the effect of ATP on CLCT. We have also con- 
structed double mutants and have studied the effects of 
ATP and CTP or ATP and GTP on the wild-type (WT) 
and mutant channels. The results from ATP and CTP or 
ATP and GTP modulations of these channels suggest 
that the residue at position 634 very likely interacts 
directly with the base moiety of the nucleotide. Further- 
more, individual mutations also show nonadditive effects 
on the ATP modulation of CLCT, a sign for an energy 
coupling between these residues in interacting with the 
bound ATP molecule. 

MATERIALS AND METHODS 

Construction of the homology model of the C-terminal 
cytoplasmic domain of CLC-1 

The sequence of human CLC-1 (Steinmeyer et al., 1991) from 
L598 to T874 was used to construct the homology model, corre- 
sponding to the sequence of CLC-5 from H578 to N737. Two 
segments in the CLC-1 C terminus, one from L598 to R672 (cor- 
responding to H578 to 1641 in CLC-5) and the other from P825 
to T874 (corresponding to P687 to A736 of CLC-5), contain the 
sequence for the CBS1 and CBS2, respectively. However, the seg- 
ment from residue 671 to residue 824 of CLC-1 is missing in 
CLC-5. It is unknown whether this segment is disordered (i.e., 
inuinsically unstructured) or not. Excluding this missing segment, 
the sequence alignment using Vector NTI software (Invitrogen) 
showed a sequence identity of 19.7% and a similarity of 31.3% be- 
tween CLC-1 and CLC-5. Homology modeling was performed 
with a bound ATP molecule using the MODELLER program 
(Marti'-Renom et al., 2000). The structures of the two segments 
containing the sequence of CBS domains of CLC-1 (Fig. 1, red 
and yellow) were determined using the crystal structure of CLC-5 
as a template (Protein Data Bank accession no. 2J9L) with energy 
minimization using the MODELLER program. It can be seen that 
the missing segment between the two CBS domains (Fig. 1, dotted 
line) is remote from the bound ATP and thereby has little 
effect on the structures of the two CBS domains forming the ATP- 
binding site (Fig. 1). 



Mutagenesis and channel expression 

The ATP modulation of CLC-1 is controlled by oxidation and re- 
duction (Zhang et al., 2008). We therefore chose Xenopus laevis 
oocytes as the channel expression system, and relatively large 
electrodes (with a tip diameter of ~7-10 um) were used to excise 
large membrane patches because the ATP effect had previously 
been shown to be more stable using these experimental approaches 
compared with those using small patches excised from mamma- 
lian cells (Zhang et al., 2008) . To further avoid decline of the cur- 
rent as a result of potential oxidation, all electrophysiological 
experiments were conducted in bath solutions containing 50-100 uM 
of reducing reagents (dithiothreitol or |3-mercaptoethanol) . 

All mutants and the WT human CLC-1 Cl~ channel were con- 
structed in the pTLN vector. Messenger RNAs were synthesized 
using the SP6 mMessage mMachine kit (Invitrogen). Harvesting 
Xenopus oocytes and the injection of messenger RNAs followed 
routine procedures as described previously (Chen, 1998; Li et al., 
2005) . Excised inside-out patch recordings were performed 3-5 d 
after RNA injections. The point mutants were constructed using 
the QuickChange II Site-Directed Mutagenesis kit (Agilent Tech- 
nologies). The mutated region was confirmed by commercial 
DNA sequencing services. 

Electrophysiological recordings 

Electrophysiological experiments were conducted using the Axo- 
patch 200B amplifier and the Digidata 1 320 A/D board controlled 
by pClamp8 software (Axon Instruments/Molecular Devices). 
Signal digitization and filtering (by the digital filter in the 
pClamp8 program) were 10 kHz and 5 kHz, respectively. The tip 
diameter of the recording electrodes was relatively large, and 
these electrodes normally had a resistance of 0.4-0.6 Mft when 
filled with a pipette (extracellular) solution containing 120 mM 
iV-methyl-D-glucamine-Cl, 1 mM MgCl 2 , 10 mM HEPES, and 1 mM 
EGTA, pH 7.4. Because the ATP effect was not prominent at a 
neutral pH, we examined the effect of ATP at pH 6.8. The bath 
(intracellular) solutions had the same ionic components, with 
pH being adjusted to 6.8 after the desired concentrations of ATP 
were added. Mg 2+ -ATP was purchased from Sigma-Aldrich. A stock 
solution of 1 00 mM was made in distilled water and was stored at 
— 20°C. Working solutions of ATP were made on the same day of 
the experiments. 

For any given patch, the membrane potential of the excised 
patch was held at 0 mV. The voltage protocol used to construct 
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Figure 1. Homology modeling of the C-terminal 
domain of CLC-1. (A) Homology model of the CBS 
domain of CLC-1 shown in ribbon presentation. The 
two segments of CLC-1 sequence corresponding to 
CLC-5 are colored in red and yellow, whereas the 
segment that is missing in CLC-5 is illustrated with a 
dotted black line. ATP is shown in stick presentation 
(in CPK color). Key amino acids examined in this 
study are colored in yellow (V613), orange (V634), 
blue (V860), and purple (E865). (B) Sequence align- 
ment of CLC-0, CLC-1, CLC-5, and AtCLC-a for the 
region encompassing the four residues shown in A. 
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the voltage dependence of the common gate open probability 
(P 0 C -V curve) has been described previously (Tseng et al., 2007; 
Zhang et al., 2008) . This voltage protocol starts with a test voltage 
varying from 120 to —140 mV (in -20-mV steps) for 400 ms fol- 
lowed by a brief depolarizing voltage to 170 mV for 400 us and 
then finally a tail voltage at —100 mV for 300 ms. Channel cur- 
rents were recorded in a given patch both in the absence and in 
the presence of ATP. The current in response to the tail voltage 
was fitted to a double-exponential function, and the initial value 
of the tail current was determined by extrapolating the fitted 
curve to the start of the tail voltage. The initial tail current of each 
trace was then normalized to the initial tail current obtained after 
the 120-mV test voltage in the absence of ATP. Because the fast 
gate open probability is fully activated by the brief voltage step at 
170 mV, the normalized initial tail current represents the open 
probability of the common gate (P 0 C ) at the preceding test voltage 
(Accardi and Pusch, 2000; Duffield et al., 2003; Bennetts et al., 
2005) . The initial tail current in the presence of ATP was also nor- 
malized to the maximal tail current (with 120-mV test voltage) in 
the absence of ATP. The normalized value (or P 0 C ) was then plot- 
ted against the test voltage, V. The half-activating voltage (Vi/ 2 ) of 
the .P 0 C -V curve was determined by fitting the data points to a 
Boltzmann equation: 

^o C = ^ mi „+ (l"^n)/{l + exp[zF(V-V 1/2 )/RT]}. 

The difference in Vi/ 2 (i.e., AVi/ 2 ) between the absence and pres- 
ence of ATP was then determined as the qualitative response to a 
certain concentration of ATP. 

Mutant cycle analyses 

Data analyses and presentations were performed using the combi- 
nation of pClamp8 and Origin software (Origin Laboratory). One 
way of evaluating the apparent affinity of ATP was to construct the 
dose-dependent curve by plotting AVi/ 2 as a function of the ATP 
concentration. The apparent K l/2 value was then determined ac- 
cording to the Michaelis-Menten equation 

AV 1/2 = AV 1/2 , max / (l + K 1/2 / [ATP]) , ( 1 ) 

where AVi/ 2imax is the maximal shift of the Vi/ 2 at saturated ATP 
concentrations, whereas Kj/ 2 is the apparent half-effective concen- 
tration in shifting the P 0 C -V curve. We used the aforementioned 
analysis for a qualitative comparison of the ATP effects. However, 
the apparent Ki/ 2 derived from the shift of V1/2 does not truly re- 
flect the half-point of ligand occupancy. To estimate the half-point 
of occupancy (or binding probability) of ATP, we used a four-state 
allosteric model to describe the opening and closing of the chan- 
nel with or without ATP binding (Zifarelli and Pusch, 2009) : 



M l K ° 

CATP-^-OATP 

(SCHEME 1) 

In Scheme I, C and O represent the closed and open states of the 
common gate, and a and c represent the gating equilibrium con- 
stants of the ATP-unbound and the ATP-bound channels. Iv and 
Ko are the binding dissociation constants of the closed-state and 
open-state channels for ATP, respectively. To determine K„ and K,., 
the values of P 0 C at various ATP concentrations, namely P 0 C (ATP), 



were normalized to the P 0 C in the absence of ATP. The normalized 
values were then plotted against ATP concentrations (see Figs. 6 
and 7) and were fit into Eq. 2 (Zifarelli and Pusch, 2009): 

P„ c (ATP) /P: (control) = 

(1 + [ATP]/K 0 )/(1 + [ATP]/K C ). 

For mutant cycle analyses, the interaction between the ligand 
and the binding site can be altered by the perturbation on the 
ligand, the binding site, or both. The free energy of the interac- 
tion, which is determined from K„ or Iv (AG 0 = RTlnK,, and AG C = 
RTlnfv), can therefore be perturbed by manipulating the ligand 
and the receptor individually or simultaneously, thus forming a 
mutant cycle as shown in Fig. 6. If the two parts changed are inde- 
pendent of each other, the free energy perturbation caused by 
the double manipulation is equal to the sum of the free energy 
perturbations from individual manipulations. Namely, 

AG(X 2 ,Y 2 )-AG(X 1 ,Y 1 )= (JJ) 
[AG(X 1 ,Y 2 )-AG(X 1 ,Y 1 )] + [AG(X 2 ,Y 1 )-AG(X 1 ,Y 1 )] ! 

where AG(XY) is the free energy of binding of the ligand Yto the 
channel X. A rearrangement of Eq. 3 leads to 

AG(X 1 ,Y 1 ) + AG(X 2 ,Y 2 ) = AG(X 1 ,Y 2 ) + AG(X 2 ,Y 1 ) 
or 

RTln[K(X 1 ,Y 1 ) x K(X 2 ,Y 2 )] = 
RTln[K(X 1 ,Y 2 ) x K(X 2> Y 1 )], 

where K(X,Y) is K,, or Kc from the experiment in channel X using 
ligand Y. The coupling coefficient, Cl, for the two perturbations in 
this mutant cycle is defined as 

Q = [K(X 1 ,Y 1 ) x k(x 2 ,Y 2 )]/[k(x 1 ,Y 2 ) XK(X 2 ,Y 1 )]. ( 4 ) 

When 0 is close to unity, there is little coupling energy between 
these two manipulations because the coupling energy of the pair 
of manipulations is given by RTlnfl => 0. On the other hand, when 
0 is significantly different from 1 , the mutant cycle analysis sug- 
gests that the two manipulations are not independent. Previous 
studies have used mutant cycle analyses to identify the interacting 
residues located in the potassium channel pore and on the pore- 
blocking toxin, respectively (Hidalgo and MacKinnon, 1995; 
Naranjo and Miller, 1996; Ranganathan et al., 1996). 

The mutant cycle approach can be modified for examining 
whether the free energy perturbations of two mutations within the 
binding site of ATP are additive. Examples are shown in Fig. 7. 
Again, an Cl value significantly deviating from unity suggests that 
these two mutations may be coupled energetically, and thus the 
two residues may interact with the bound ligand simultaneously. 

Online supplemental material 

Fig. SI shows a comparison of the two homology models of CLGl's 
Gterminal domain based on the high resolution structures of CLG5 
and CmCLC. Fig. S2 shows the values of 0. calculated from the P c c 
values at other voltages. Online supplemental material is available at 
http://wwwjgp.org/cgi/content/full/jgp.201010495/DCl. 

RESULTS 

The C-terminal cytoplasmic regions of all eukaryote 
CLC molecules contain two tandem CBS domains that 
form a putative ATP-binding site (Scott et al., 2004; 
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Wellhauser et al., 2006). This predicted ATP-binding 
site has been direcdy observed in the high resolution 
structure of the C terminus of CLC-5 (Meyer et al., 2007) . 
So far, functional consequences of ATP modulation 
have been shown in CLC-1 as well as in CLC-2, CLC-5, 
andAtCLG-a (Niemeyer etal., 2004; De Angeli etal., 2009; 
Zifarelli and Pusch, 2009) . In this study, we use CLC-1 to 
study the interaction of ATP with CBS domains. 

We first constructed a homology model of CLC-1 's 
C-terminal region based on the high resolution, ATP- 
bound structure of the C-terminal domain of CLC-5 
(Fig. 1 A) . An independent model was built later based 



Figure 2. Effects of alanine mutations for 
four residues on the ATP modulation of 
CLC-1. (left and middle) Recording traces 
of WT CLC-1 and the four mutants. The 
recordings in the absence (left) and in 
the presence (middle) of 3 mM ATP were 
obtained from the same patch (pH; 6.8). 
(right) Steady-state P Q C -V curves of WT 
CLC-1 and the four alanine mutants shown 
in left and middle panels. Black squares, 
red circles, and green triangles represent 
the data in control (before ATP) , 3 mM 
m ' mV ' ATP, and washout, respectively. 

on the recently solved structure of CmCLC (Fig. SI A) . 
The homology model based on the CLC-5 structure 
suggested several residues that might interact with the 
bound ATP molecule, including V613, V634, V860, and 
E865. Other nearby residues, such as T636, P638, H847, 
and L848, appear to line the ATP-binding pocket as 
well. However, because the latter four residues in CLC-1 
have been examined previously (see Fig. SI B for the 
residues in CLC-1 whose mutations affect the ATP 
effect) , we focused on the former four residues in this 
study. Fig. 1 B shows the alignment of the amino acid 
sequence from CLC-0, CLC-1, CLC-5, and AtCLC-a in 
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TABLE I 

V1/2 and P min of the Pf-V curves of various mutants 
Channels/mutants V1/2 Pmin n 

WTCLC-1 -45.9±2.8 a 0.35 ± 0.02 a 15 

V634A -59.3 ±0.8 0.64 + 0.02 7 

V634D -34.0 ±3.9 0.28 ± 0.02 7 

V634F -56.8 ±2.0 0.40 ±0.01 9 

V634Y -72.4 ± 1.8 0.50 ± 0.02 7 

V634C -44.4 ±2.5 0.30 ± 0.01 7 

V634S -58.3 ±0.9 0.42 ± 0.02 4 

V634G -47.6 ±2.8 0.32 ± 0.01 7 

V634W -54.7 ±2.6 0.45 ± 0.02 4 

E865A -68.5 ±2.6 0.52 ± 0.01 4 

E865D -36.2 ±2.5 0.24 ± 0.01 11 

E865S -43.7 ±8.2 0.32 ± 0.01 3 

V860A -74.6 ± 1.3 0.67 ± 0.02 3 

V860S -54.9 ±7.6 0.26 ± 0.02 4 

V860Y -37.5 ±8.1 0.28 ± 0.01 4 

V613A -49.1 ± 3.4 0.37 ± 0.01 7 

V613F -43.9 ±8.7 0.28 ± 0.02 3 

V613G 88.2 ±4.7 0.15 ±0.01 5 

a Error from curve fitting. 

regions containing the target residues. Among these 
four CLC molecules, functional ATP modulation effects 
have been observed in the latter three CLCs: ATP inhib- 
its the activity of CLC-1 and AtCLC-a, whereas it po- 
tentiates the activity of CLC-5. The homology model 
suggested that the adenine ring of ATP appears to be 
surrounded by the side chains of V634, V860, and V613, 
whereas E865 may interact with the ribose of ATP. 
Among these four positions, two corresponding residues 



in CLC-5 were examined previously by others (Meyer 
et al, 2007; Zifarelli and Pusch, 2009). The ATP bind- 
ing was shown to be eliminated by the Y617A and D727A 
mutations of CLC-5 (corresponding to V634A and E865A 
of CLC-1). Mutating two equivalent residues, H620 
and D753 (corresponding to V634 and E865 of CLC-1) , 
also diminished the ATP inhibitory effect in AtCLC-a 
(De Angeli et al., 2009). In addition, T596 and 1722 of 
CLC-5 (corresponding to V613 and V860 of CLC-1) 
have also been suggested to be involved in ATP binding 
to CLC-5 (Meyer et al., 2007). Thus, V613, V634, V860, 
and E865 of CLC-1 are likely to interact with the bound 
ATP molecule. 

We started our experiments by examining the ATP 
modulation effects on the mutants in which the WT res- 
idues at these four positions were replaced by alanine. 
Fig. 2 shows the effects of ATP on WT CLC-1 and those 
of the four alanine mutants at pEL 6.8. The recording 
traces were elicited by the voltage protocol used for 
constructing the common gate P 0 C -V curve (see Materi- 
als and methods) . The comparison of the recording 
traces both in the absence (Fig. 2, left) and in the pres- 
ence of 3 mM ATP (Fig. 2, middle) shows a robust inhi- 
bition of WT CLC-1 by ATP (Fig. 2, top), but the effect 
of ATP on these alanine mutants is either weak or com- 
pletely lost. The P 0 C -V curves shown in Fig. 2 (right) 
indicate that the V634A, V860A, or E865A mutation 
eliminates the ATP modulation effect and the V613 
mutation significandy reduces the ATP effect. 

As an initial effort to explore the interaction of ATP 
with the potential binding site more systemically, we re- 
placed the WT residues with multiple amino acids at 
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Figure 3. Changes of the V1/2 of the 
P 0 C -V curves of various mutants in re- 
sponse to 3 mM ATP at pH; 6.8. The 
P 0 C -V curves of each mutant in the ab- 
sence and in the presence of 3 mM 
ATP were constructed from the same 
patch, and the Vi/ 2 of the P a c -V curves 
were obtained by fitting the data points 
to a Boltzmann equation. The differ- 
ence in V1/9 between control and 3 mM 
ATP (AV1/2) was nrs t determined from 
each patch and then averaged (n = 14 
for WT CLC-1 and n = 4—7 for various 
mutants). The label n.c. indicates that 
the mutant did not show functional 
current in at least two runs of expres- 
sion in Xenopus oocytes. Error bars rep- 
resent SEM. 
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Figure 4. The V634F mutation increases the 
apparent affinity for the ATP modulation. 
Representative recording traces were ob- 
tained from an excised patch of the oocyte ex- 
pressing the V634F mutant in control (0 ATP) 
and in various concentrations of ATP as indi- 
cated. The averaged P 0 C -V curves of the V634F 
mutant in various ATP concentrations are 
shown in the bottom right panel. Error bars 
represent SEM. n = 7. 



these four positions. Mutating V613, V860, and E865 
affects the functional expression; at each site, four out 
of seven mutants did not generate functional current in 
Xenopus oocytes. However, mutations at the V634 posi- 
tion are more tolerable; all mutants made in this study 
at this position showed functional current. Table I sum- 
marizes the minimal open probability (P m i n ) and the 
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Figure 5. Shift of the V1/2 of the P 0 C -V curve (AV1/2) m various 
mutants of V634 as a function of the ATP concentration. Data 
points were fitted to Eq. 1 to determine the apparent K 1/2 and 
Vi/ 2 ,max. (A) Mutants of V634 in which the apparent ATP affinity 
is reduced. For comparison, the ATP dose-response curve of WT 
CLC-1 is also shown. (B) ATP dose-dependent curves of WT CLC-1, 
V634F, V634Y, and V634W. Error bars represent SEM. n = 4-8. 



half-activation voltage (V 1/2 ) of the common gate P Q C -V 
curve of these functional mutants. We then tested the 
shift of the P Q C -V curve in the functional mutants by 
3 mM ATP at pH ; 6.8. Fig. 3 shows the change of the V 1/2 
(AV1/2) of the P Q C -V curve in response to 3 mM of intra- 
cellular ATP (pH ; 6.8) . For mutations at positions V613, 
V860, and E865, all functional mutants show a sub- 
stantial reduction of AV1/2 in response to 3 mM ATP. 
However, the effects of mutations at the V634 position 
are diverse. The V634A, V634D, V634G, and V634S muta- 
tions suppress the ATP effects, whereas the ATP modu- 
lation of V634C is similar to that of the WT channel. 
Interestingly, some mutants, for example V634F, show a 
higher sensitivity to the ATP modulation. Although the 
AV1/2 of V634F generated by 3 mM ATP appears to be 
similar to that of WT CLC-1 (Fig. 3) , the apparent affin- 
ity of ATP in shifting the P Q C -V curve of this mutant is 
higher than that of WT CLC-1 — the shift nearly reached 
a maximum at a concentration of 0.1 mM (Fig. 4). 

To examine the diverse functional consequences of mu- 
tating residue V634, we studied the ATP dose-dependent 
effect in various V634 mutants. Fig. 5 A shows the dose- 
response curves of those mutants with less ATP sensitiv- 
ity, whereas Fig. 5 B depicts the results of mutants with a 
higher ATP sensitivity. The K 1/2 and the maximal value 
of the V1/2 shift (AVi /2jimx ) are shown in Table II. The 
ATP modulations of V634C and V634S mutants are 
roughly similar to that of the WT channel: the AV 1/2j max 
in both cases is similar to that of the WT channel, 
whereas the apparent affinity is only slightly reduced in 
the V634S mutant. However, the V634D and V634G mu- 
tants show a more severe disruption in ATP regulation: 
the apparent ATP affinity and/or the maximal shift of 
the V1/2 is reduced in these two mutants. Finally, the 
V634A mutation completely eliminates the ATP modu- 
lation effect; even at an ATP concentration of 10 mM, 
no significant shift of V 1/2 is observed. However, placing 
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TABLE II 



Ki/2 and A Vi/2 1 mafo r the ATP dose— response curves ofV634 mutants 



Channels/ mutants 






1} 




mM 


mV 




WT CLC-1 


0.78 ± 0.08 


76.9 ±7.4 


8 


V634F 


0.026 ± 0.002 


61.6 ±4.9 


7 


V634Y 


0.050 ± 0.004 


46.6 ± 3.3 


5 


V634W 


0.11 ±0.01 


90.7 ±9.4 


4 


V634C 


0.81 ±0.06 


70.8 ± 3.3 


5 


V634S 


1.8 ±0.3 


72.9 ± 8.5 


4 


V634D 


7.4 ± 1.5 


18.5 ±3.1 


4 


V634G 


13.3 ±4.0 


77.8 ±8.1 


4 


E865D 


0.40 ± 0.09 


14.2 ±2.4 


6 


V613A 


2.3 ± 0.3 


26.3 ± 3.4 


4 


E865D-V634F 


0.089 ± 0.05 


8.7 ±2.4 


4 


V613A-V634F 


2.3 ± 0.35 


62.6 ± 5.4 


7 



K 1/2 determined from Eq. 1, using AV 1/2 as the response to ATP inhibition. 
Means ± SEM are given. 

aromatic amino acids at position 634 appears to en- 
hance the ATP modulation. The apparent Ki/ 2 of the 
ATP effects are 0.026, 0.050, and 0.11 mM for V634F, 
V634Y, and V634W (Table II) , corresponding to an in- 
crease of the apparent ATP affinity by °^30-, ^16-, and 
^^-fold, respectively. 

Because only those mutants with an aromatic side 
chain at the V634 position showed a high ATP sensitivity 
(namely V634F, V634Y, and V634W) , we suspected that 
the ring-ring interaction may contribute to the higher 
apparent ATP affinity in these mutants. If the adenine 
of ATP indeed direcdy interacts with the side chain of 
residue at position 634, we expect that a mutant cycle 
analysis should reveal that the free energy perturbations 
caused by mutating residue 634 and by altering the base 
of the nucleotide should couple energetically. 

However, the estimate of the apparent affinity of ATP 
inhibition based on the change of V1/2 could be prob- 
lematic because AV1/2 m ay not truly reflect the change 
of ligand occupancy. For this reason, we model the ATP 
inhibitory effect using a four-state allosteric scheme 



(Scheme I) to evaluate the affinities of WT CLC-1 and 
V634F mutants for ATP and other nucleotides. Fig. 6 
gives an example of mutant cycle analysis in which the 
normalized P Q C at —40 mV in different conditions are 
plotted against ligand concentrations. Fig. 6 A shows 
that the free energy changes caused by individual per- 
turbations are not additive, with large fl values (see 
Materials and methods) of 60 and 113 calculated from the 
values of Kc and Ko, respectively. We have also performed 
similar analysis for WT CLC-1 and V634F using ATP and 
GTP as ligands (Table III) , with a calculated fl of ^7.6-10. 
However, when deoxy-ATP (dATP) was used to examine 
the interaction of the base with the side chain of resi- 
due 634, O values close to 1 were obtained (Fig. 6 B 
and Table III). Analyses of P Q C at various voltages re- 
vealed similar fi values in the voltage range between 
— 140 and —40 mV (Fig. S2 A). These results again are 
consistent with the idea that the base of the bound nu- 
cleotide directly interacts with the side chain of residue 
634 of CLC-1. 

The homology model of CLC-1 also suggests that 
V613, V860, and E865 may interact with ATP. We thus 
examined the perturbations of the binding free energy 
of ATP to see whether these changes induced by indi- 
vidual mutations are additive. Unfortunately, construc- 
tion of the ATP dose-dependent curve for V860 mutants 
is difficult because of the very small shift of the P a -V 
curve caused by ATP. We therefore focused on V613 
and E865 and studied the synergistic effects between 
V634F and V613A and between V634F and E865D mu- 
tations because these mutations still retain some degree 
of ATP inhibition and at the same time alter the appar- 
ent ATP affinity. The concentration-dependent curves 
of the ATP-induced change of P Q C for the single and 
double mutants were constructed, and the ATP affinity 
(Kq and Kc of Scheme I) for each mutant was compared 
with that of WT CLC-1. Fig. 7 A shows a double mutant 
cycle analysis that examines the interaction of the V634F 
and V613A mutations, whereas Fig. 7 B shows the analy- 
sis of the interaction between V634F and E865D. 



TABLE III 

Apparent affinities of ATP, CTP, GTP, and dATP for WT CLC-1 and V634F 



Channels/ mutants Ko or Kc 





ATP 


CTP 


GTP 


dATP 




mM 


mM 


mM 


mM 












WT CLC-1 


0.32 ±0.04 (8) 


1.1. ±0.2 (6) 


9.9 ± 1.0 (5) 


0.85 ±0.14 (4) 


V634F 


0.017 ±0.005 (7) 


3.5 ± 0.4 (8) 


4.0 ± 0.6 (6) 


0.049 ±0.013 (6) 






fi = 60 


a = 7.6 


n = 1.09 


K„ 










WT CLC-1 


0.70 ±0.10 (8) 


1.3. ±0.3 (6) 


15.1 ± 1.9 (5) 


1.3 ±0.2 (4) 


V634F 


0.039 ±0.013 (7) 


8.2 ±1.2 (8) 


8.5 ± 1.6 (6) 


0.071 ±0.021 (6) 






£1 = 113 


a = 10 


0 = 0.98 


Calculated from the values of P a c at 


— 40 mV. Means ± SEM are given. 


The numbers inside parenfh 


sses indicate the 


number of patches. 
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Figure 6. Interaction of nucleotides with residue V634 in the ATP-binding site explored by double mutant cycle analyses. (A) Mutant 
cycle analysis of the interaction of the nucleotide with residue V634. The ATP and CTP dose-response curves of WT CLC-1 and the 
V634F mutant were first constructed using the parameter P 0 c (ATP)/.P 0 c (control) as the dependent parameter. The values of Ko and Kc 
(Table III) of ATP/CTP modulations were determined by fitting the data points to Eq. 2, and the Kc values were used to calculate the 
free energy differences of single or double manipulations (bar graph) according to Eq. 3 and the coupling constant, CI, according to 
Eq. 4. n = 6-8. (B) Mutant cycle analysis of the nucleotide interaction with V634 using ATP and dATP. The ligand dose-response curves 
of WT CLC-1 and V634F were constructed using P D C (ATP) /P a c (control) as the dependent parameter. Ko and Kc were determined (Table III), 
and the fl value and the free energy differences were calculated accordingly, n = 4—8. 



Although the value of fl for the former pair of muta- 
tions appears to be large, whereas that of the latter pair 
of mutations is smaller (see Fig. S2 B for fl values at dif- 
ferent voltages) , the results in both cases suggest that 
the mutational effects of V634F and those of V613A and 
E865D on the ATP regulation may couple together, a 
result consistent with the idea that these residues inter- 
act with ATP simultaneously. 

DISCUSSION 

CLC-1 CI channels contribute a major conductance to 
the sarcolemmal membrane of skeletal muscles and, 
thus, are critical for the membrane excitability of skele- 
tal muscles. The importance of CLC-1 in maintaining 
normal muscle functions can be seen from the myoto- 
nia disease resulting from defects of CLC-1 (Koch et al., 
1992). Likely, the CLC-1 conductance facilitates repo- 
larization of the membrane of skeletal muscles after firing 



an action potential (Jentsch et al., 2005). In addition, 
recent studies suggested that the reduction of the CI 
conductance of the sarcolemma by a low intracellular 
pH (pH;) can help overcome muscle fatigue (Pedersen 
et al., 2004, 2005). Experiments performed several de- 
cades ago showed that acidification of muscle cells can 
reduce the muscle CI conductance (Hutter and Warner, 
1967a,b) . However, various experiments on the cloned 
CLC-1 channel showed that reducing pH ; increases the 
open probability of CLC-1 (Rychkov et al., 1996; Saviane 
et al., 1999; Accardi and Pusch, 2000). On the other 
hand, we and others recently found that if ATP is in- 
cluded in the intracellular solution, reducing pH; de- 
creases the CLC-1 activity (Bennetts et al., 2007; Tseng 
et al., 2007). Thus, binding of ATP to the intracellular 
side of CLC-1 appears to be a prerequisite for reducing 
the CLC-1 conductance by a low pH ; . 

The primary sequences of CLC channels/ transporters 
have suggested the presence of a potential ATP-binding 
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Figure 7. Nonadditive effects of mutations on the ATP modulations of CLC-1. The parameter P 0 C (ATP) / P 0 C (control) was used as the 
response to ATP modulation. (A) Mutant cycle analysis of the effects of V634F and V613A mutations on the ATP modulation. ATP 
dose-response curves of single point mutants and the double mutants were constructed to determine the values of Ko and The fitted 
K,- values were used to calculate the free energy differences of single or double mutants (bar graph) according to Eq. 3 and the coupling 
constant, fl, according to Eq. 4. In the context of Eq. 2, the V634F mutation corresponds to the change from Xi to X 2 , whereas the V613A 
mutation represents the change from Yi to Y 2 . (B) Double mutant cycle analysis for V634F/E865D mutations. The values of Cl at different 
voltages are shown in Fig. S2. Error bars represent SEM. n = 4—8. 



site in the CBS domain located at the C terminus of 
CLC molecules. Structural and biochemical experiments 
on the C-terminal domain of CLC-5 demonstrated that 
ATP indeed binds to the predicted CBS domain of CLC-5 
(Meyer et al., 2007). Does ATP also bind to the CBS 
domain of CLC-1 and cause the inhibition of the CLC-1 
current? Recent studies from different laboratories pro- 
vided somewhat inconsistent results. Our results (Tseng 
et al., 2007) and those from Bennetts et al. (2007) 
showed that ATP may modulate the common gating of 
CLC-1. On the other hand, Zifarelli and Pusch (2008) 
failed to show an ATP effect on CLC-1 and thus sug- 
gested that the ATP modulation of CLC-1 may be 
achieved via an accessory protein. Although our further 
experiments showed that the ATP regulation of CLC-1 
is modulated by oxidation and reduction (Zhang et al., 
2008) and thus suggested that the variation of the oxida- 
tion reduction level in various experimental conditions 
could be one of the underlying reasons for inconsistent 
experimental findings, the evidence of a direct ATP 



binding to CLC-1 has still been lacking. With the 
studies outlined here, we have direcdy examined the 
interaction of ATP with the amino acid residues lining 
the potential ATP-binding pocket in the CBS domains 
of CLC-1. 

To identify the potential residues in CLC-1 that may 
interact with the ATP molecule, we constructed homol- 
ogy models of CLC-1 based on the C-terminal structures 
of CLC-5 and CmCLC. Previously, two homology mod- 
els have been reported for the C-terminal domain of 
CLC-1. One of them (Bennetts et al., 2005) was con- 
structed based on the crystal structure of inosine-mono- 
phosphate dehydrogenase (Zhang et al., 1999) , which is 
not very reliable because of the low sequence identity 
between CLC-1 and inosine-monophosphate dehydro- 
genase (Bennetts et al., 2007). The other model rebuilt 
by Bennetts et al. (2007) was based on the crystal struc- 
ture of the C-terminal domain of CLC-0 (Protein Data 
Bank accession no. 2D4Z), which is an apo structure 
with no bound ATP molecule. The alignment between 
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the C-terminal structures of CLC-0 and CLC-5 shows a 
somewhat different orientation of the two CBS domains 
(Meyer et al., 2007). Similarly, the two CLC-1 C-terminal 
models we built based on CLC-5 and CmCLC structures 
also show quite a significant difference in the relative 
positions of the two CBS domains (and thus the resi- 
dues interacting with ATP; Fig. SI) . It cannot be excluded 
at this moment that part of the difference among vari- 
ous C-terminal structures results from an ATP-induced 
conformational change. Because the structure of the 
CLC-5 's C-terminal domain contains ATP, we used 
the homology model based on this structure to guide 
our experiments. 

Although CLC-1 and CLC-5 only show ^20% sequence 
identity, the homology model of the potential ATP- 
binding site in CLC-1 is quite similar to that in the 
C-terminal structure of CLC-5. For example, the resi- 
dues V634 and E865 that correspond to Y617 and D727 
of CLC-5 in sequence alignment are also located in the 
ATP-binding sites, and their side chains appear to inter- 
act with the adenine ring and the ribose of the bound 
ATP, respectively (Fig. 1 ) , just like in the structure of the 
CLC-5 C-terminal domain. Mutating both V634 and 
E865 into alanine eliminates the ATP effect, a result 
consistent with biochemical data from CLC-5. In addi- 
tion, the homology model also suggests that two other 
residues, V613 and V860, may also interact with ATP. 
Mutating these four residues in the potential ATP- 
binding pocket resulted in various consequences. 
Overall, the mutations at the V634 position are relatively 
tolerable for the functional expression of the channel. 
All mutations here resulted in functional mutant chan- 
nels. For the mutations at positions V613, V860, and 
E865, four out of seven mutations at each position did 
not produce functional current. We have not character- 
ized whether the mutations at these three positions 
resulted in nonfunctional channels or whether the mu- 
tations simply caused a problem of membrane traffick- 
ing. For those mutants that can generate functional 
current, the ATP regulation of CLC-1 was altered. 

Overall, the majority of the mutants cause a suppres- 
sion of the ATP effect; either the apparent affinity (K 1/2 ) 
or the efficacy (the maximal AV1/2) of the ATP modula- 
tion is reduced, or both. However, introducing an aro- 
matic residue at position 634 generates mutant channels 
with a higher apparent ATP affinity. It is worth pointing 
out that WT CLC-1 contains a valine residue at position 
634, whereas the corresponding residue in CLC-5 is 
tyrosine. The apparent K 1/2 of the ATP in WT CLC-1 at 
pH ; 6.8 is higher than that in CLC-5 (^90 uM) observed 
by Meyer et al. (2007) based on a binding assay, al- 
though a patch-clamp study by Zifarelli and Pusch (2009) 
reported that the ATP affinities of CLC-5 are 0.43 mM 
and 0.89 mM for the active and inactive states of this 
transporter, respectively. Nonetheless, both studies showed 
that the corresponding residue (Y617) in CLC-5 is 



crucial for the ATP modulation. Because placing aro- 
matic amino acids at position 634 increased the appar- 
ent ATP affinity, we hypothesized that the aromatic ring 
of these amino acids may directly interact with the ade- 
nine of ATP. To further explore this possibility, we stud- 
ied the apparent affinities of various nucleotides on 
WT CLC-1 and on the V634F mutant and showed that the 
O values from mutant cycle analyses are significandy de- 
viated from unity for ATP/CTP and ATP/GTP pairs 
(Fig. 6 A and Table III) . However, using ligands that dif- 
fer in the ribose moiety (dATP vs. ATP) as the modula- 
tion ligands showed that the O value is close to unity 
(Fig. 6 B and Table III) . These results are consistent with 
the hypothesis that the base moiety of the bound nucle- 
otide in the CBS domain of CLC-1 interacts directly with 
residue V634. 

Thus, the amino acid residue of V634 in CLC-1, like 
the corresponding residue in CLC-5 (Y617), should be 
critical for ATP binding to the CBS domain. In addition 
to V634, three other amino acid residues shown in this 
study also affect the ATP modulation of CLC-1 . We won- 
dered whether the mutations at these various positions 
affect the ATP binding independendy. If so, the free en- 
ergy change generated by the double mutation is ex- 
pected to be a summation of the free energy changes 
from individual mutations. For this purpose, we used 
mutant cycle analyses to examine the interaction between 
the mutation V634F and two other mutations, V613A 
and E865D. We chose these two latter mutations be- 
cause the mutants still retain some ATP inhibitions for 
constructing a reliable ATP dose-dependent curve. The 
apparent K 1/2 of the ATP modulations for the single 
point mutants and the double mutants are shown in 
Table II. Fig. 7 A shows the double mutant cycle analysis 
for V613A and V634F mutations. The results show a 
violation of the independency of the two mutations be- 
cause the fi value from this mutant cycle analysis clearly 
deviates from unity, indicating that the free energy 
change of the double mutant V613A/V634F (compared 
with the WT channel) is not the sum of the free energy 
changes from the V613A and V634F mutations. In com- 
parison, experiments with V634F and E865D mutants 
appear to have a smaller fl value of only ^6.7 (Fig. 7 B) . 
We have no definitive interpretation of the large cou- 
pling energy between V613A and V634F mutations and 
the small coupling energy between V634F and E865D 
except speculating that V6 13 and V634 interact with the 
adenine ring of ATP, whereas V634 and E865 may be 
separated in farther distance. Nonetheless, the nonad- 
ditive effects of individual mutations and the strong en- 
ergetic coupling in ATP/CTP and ATP/ GTP bindings 
provide evidence against the idea that the ATP modula- 
tion on CLC-1 is mediated by an accessory protein. 

Previous studies by others and our own experiments 
showed that the inhibition of the CLC-1 common gat- 
ing by ATP can also be achieved by ADP and AMP 
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(Bennetts et al., 2005; Zhang et al., 2008). It was shown 
that there is no discrimination in binding ATP, ADP, 
and AMP to the CBS domains of CLC-5 (Meyer et al., 
2007) . This similarity seems to argue for a similar overall 
architecture of the two binding sites, although the 
sequence homology of the CBS domains of CLC-1 
and CLC-5 is quite low (~30%). However, binding of 
nucleotides/nucleosides to CLC-5 potentiates but does 
not inhibit the transport activity of CLC-5 (Zifarelli and 
Pusch, 2009) . In AtCLC-a, however, it was found that 
ATP has an inhibitory effect, whereas ADP and AMP do 
not have any direct effect except that AMP is able to 
compete with ATP for binding (De Angeli et al., 2009) . 
The molecular basis underlying these detailed functional 
variations among the nucleotide/ nucleoside effects of 
these three CLC molecules remains unclear. One possible 
explanation is that binding of nucleotides/nucleosides to 
the CBS domain may induce different allosteric effects 
in different CLC molecules. Thus, the ATP modulation 
mechanisms could be quite different in various CLC 
molecules (Zifarelli and Pusch, 2009). 

In summary, the experiments presented in this study 
have identified several amino acid residues that may di- 
rectly interact with the bound ATP molecule in the CBS 
domain of CLC-1. The results thus argue that ATP mod- 
ulates the common gating of CLC-1 via a direct binding 
to the C-terminal domain of CLC-1. Given the large 
shift of the P 0 C -V curve, modulation of the CLC-1 com- 
mon gating should be a useful model to study ATP inter- 
action with the CBS domain. It will also be interesting to 
examine the functional roles of the corresponding 
residues in other CBS domain-containing molecules. 
Such a comparison of the effects of the corresponding 
mutations in different molecules could provide a deeper 
understanding of the molecular interaction between 
nucleotide and the CBS domain. 
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